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Abstract We discuss the theoretical framework to describe quasiparticle elec- 
tric and heat currents in NIS tunnel junctions in the dirty limit. The approach is 
based on quasiclassical Keldysh-Usadel equations. We apply this theory to dif- 
fusive NIS'S tunnel junctions. Here N and S are respectively a normal and a su- 
perconducting reservoirs, I is an insulator layer and S' is a nonequilibrium super- 
conducting lead. We calculate the quasiparticle electric and heat currents in such 
structures and consider the effect of inelastic relaxation in the S' lead. We find 
that in the absence of strong relaxation the electric current and the cooling power 
for voltages eV < A are suppressed. The value of this suppression scales with 
the diffusive transparency parameter. We ascribe this suppression to the effect of 
backtunneling of nonequilibrium quasiparticles into the normal metal. 

Keywords Electron cooling • NIS tunnel junctions 
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1 Introduction 

It is well known that in NIS (Normal metal - Insulator - Superconductor) tunnel 
junctions the flow of electric current carried by quasiparticles is accompanied by 
a heat transfer from the normal metal into the superconductor 1 1, 2]. This happens 
due to the presence of the superconducting energy gap A, which induces selective 
tunneling of high-energy quasiparticles out of the normal metal. In the tunnel- 
ing event only quasiparticles with energy E > A (compared to the Fermi level) 
can tunnel out of the normal metal. They generate the single particle current and 
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associated heat current. At lower energies E < A charge transfer occurs via mech- 
anism of Andreev reflection |3l|4[j. Andreev current Ia generates a Joule heating 
IaV that is deposited in the normal metal electrode but this effect dominates over 
single particle cooling only at very low temperatures |5]. For the temperatures 
considered in this paper Andreev Joule heating is negligible. 

It has been shown that for voltage-biased NIS tunnel junctions heat current out 
of the normal metal (also referred as "cooling power") is positive when eV < Z\, 
i.e., it cools the normal metal |6]. For eV > A the current through the junction 
increases strongly, resulting in Joule heating IV and making the heat current neg- 
ative. The cooling power is maximal near eV Ki A. 

This effect enables the refrigeration of electrons in the normal metal. A mi- 
crorefrigerator, based on an NIS tunnel junction, has been first fabricated by Na- 
hum et al. iH]. They have used a single NIS tunnel junction in order to cool a small 
normal metal strip. Later Leivo et al. |6] have noticed that the cooling power of 
an NIS junction is an even function of an applied voltage, and have fabricated a 
refrigerator with two NIS tunnel junctions arranged in a symmetric configuration 
(SINIS), which gives a reduction of the electronic temperature from 300 mK to 
about 100 mK. This significant temperature reduction gives a perspective to use 
NIS junctions for on-chip cooling of nano-sized systems like high-sensitivity de- 
tectors and quantum devices [TJ . 

To enhance the performance of the NIS refrigerator it is important to under- 
stand the role of possible factors that may facilitate or decrease the cooling effect. 
One of such factors is the inelastic relaxation of injected quasiparticles in the su- 
perconducting lead. In non-reservoir geometries the quasiparticles injected into 
the superconducting lead generate a nonequilibrium distribution. In a diffusive su- 
perconductor backscattering on impurities and subsequent backtunneling into nor- 
mal metal may considerably reduce the net heat current out of the normal metal 
electrode. 

The purpose of this work is to investigate the importance of the nonequilib- 
rium quasiparticle distribution and consider the effect of inelastic relaxation in 
the superconductor. Possible mechanisms of inelastic relaxation could be usual 
processes, such as electron-electron or electron-phonon interactions, but also the 
presence of so called "quasiparticle traps", i.e. additional normal metal electrodes 
connected to the superconductor, which remove excited nonequilibrium quasipar- 
ticles from the superconducting lead |8, 9]. In this paper we do not specify the 
mechanisms of inelastic relaxation and consider the relaxation time approxima- 
tion approach. Effect of the nonequilibrium quasiparticle injection in NIS junc- 
tions was also discussed in Ref. iflOll . where the authors proposed a phenomeno- 
logical model of quasiparticle diffusion in the superconducting lead. 

The paper is organized as follows. In the next section, we formulate the theo- 
retical model and basic equations. In sections [3] and |4] we solve the kinetic equa- 
tions and apply solutions for the calculation of the electric and heat currents, re- 
spectively. In Sec.|5]we calculate quasiparticle distribution functions in the super- 
conducting lead. Finally we summarize the results in Sec.|6l 
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Fig. 1 Geometry of the considered system. 



2 Model and basic equations 

The model of the N-I-S'-S junction under consideration is depicted in Fig. [T] and 
consists of a vokage-biased normal metal reservoir (N), an insulator layer (I), a su- 
perconducting layer (S') of thickness L and a superconducting reservoir (S) along 
the X direction. The S and S' leads are made from the identical superconducting 
material. We assume the S'-S interface to be fully transparent. We will consider 
the diffusive limit, in which the elastic scattering length £ is much smaller than 
the coherence length = ^'I'/2zi, where D is the diffusion coefficient (we as- 
sume h = ks = 1). The length L of the S' lead is assumed to be much larger than 
^0- The problem of current flow through diffusive N-I-S'-S structures with short 
S' superconductor lead was solved in Ref. [11]. 

Under the conditions described above, the calculation of the electric and heat 
currents requires solution of the one-dimensional Keldysh-Usadel equations fl3h 
(see also review Ll4i1 ) for the 4x4 matrix Keldysh-Green function G{x,E) in the 
S' lead, 

[a,E + A, G] =iDdJ, J = GdG, G^ = l, (1) 
^G^^ G'<=ff-fg\ (2) 



Here 



^ 



a, 0\ X _ M 
aj' ^-[oA 



^RA jjjg 2x2 Nambu matrix retarded and advanced Green's functions, / = 
/+ + Oyf- is the matrix distribution function (we use "check" for 4x4 and "hat" 
for 2 X 2 matrices), ^ are the Pauli matrices in the Nambu space, A = e"^~^iaxA, 
A and X are the modulus and the phase of the pair potential, and d = 3/ dx. In 
Eqs. ([Til we neglect the inelastic collision term which will be taken into account 
later. Since we are interested in small voltages eV <A (when the cooling of the 
normal metal occurs) we can neglect the effect of the suppression of the supercon- 
ducting gap due to the heating in the superconductor. 

The boundary conditions for the function G and the matrix current / at the 
left normal (x = — 0) and the ri ght superconducting {x = +0) sides of the tunnel 
barrier are given by the relation lllSll 



{aN/gN)J-o = J+o = {W/^o) [G-o,G+o] , 



(3) 



4 



where On and g^i are the normal conductivities of the N and S' leads per unit 
length, respectively. In Eq. ([3]l, the transparency parameter W is defined as 

W = R{^o)/2R = momr » r, (4) 

where R{^o) = ^o/8n is the normal resistance of the S' lead per length and R is 
the junction resistance. It has been shown in Refs. |[l^[l7l] that this quantity rather 
than the barrier transparency F plays the role of a transparency parameter in the 
theory of diffusive tunnel junctions (see also the discussion in Ref. [18]). In this 
paper, we will consider the limit W ^ 1, which corresponds to the conventional 
tunneling concept. At the right transparent S'-S interface all functions and their 
first derivatives are to be continuous. We neglect possible small resistance at this 
interface il2ll since it is much smaller than the resistance of the S' lead in the 
normal state. 

The electric and heat currents are related to the Keldysh component of the 
matrix current / respectively as lfl3l [l9 ll20l l2lll22h. 

1=-^ f TicrJ^dE, (5) 

4e Jq ' ' 

P = IV+^ [ ETrf^dE, (6) 
Jo 

and thus they can be expressed through the boundary value y+o in Eq. (O, 

I=-:^J^^TrG,[G-o,G+ofdE, (7) 

P = '^ + -^l ETr[G ^0,0+0]" dE. (8) 

For further consideration it is convenient to write the Green function in the 
following standard way, 

g = a,u + iayV. (9) 

Here we neglect the phase of the anomalous Green function since it gives correc- 
tions to the next order in diffusive barrier transparency parameter W. The func- 
tions u and V determine the spectral characteristics of the system. In particular, 
the quantity N{E) = {u^ — i/-^ /2 is the density of states (DOS) normalized to 
its value Nf in the normal state. In what follows, we will express the advanced 
Green functions through the retarded ones, (m, v)'* = — (m, v)^*, using the general 
relation^ = —a^g^^a^, and omit the superscript/?, considering retarded Green's 
functions only. 

In this paper we will calculate the single particle current, therefore we consider 
quasiparticle energy E > A from now on. We neglect the proximity effect since 
proximity corrections to the spectral functions are of the order of W . Therefore, 
we use the BCS density of states in both S' and S layers. In the left voltage-biased 
normal metal reservoir (N) we have. 



^ , E + eV\ , ^ fE-eV 
tanh — ± tanh 



27a/ / V 27) 



(10) 
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where Tn is the temperature of the normal metal reservoir. We assume the voltage 
V to be directly applied to the tunnel barrier and neglect a small electric field 
eVW) penetrating the S' superconducting lead. 

In the right superconducting reservoir (S) Green's and distribution functions 
are given by the relations 

iE,A) 

gs = o,us + iayVs, ius,vs) = ^=====, (11) 

/+5L>,^^ = tanh(A^, f_sl^=0, (12) 

where 7$ is the temperature of the S reservoir. 

In the S' layer, Green's function is given by Eq. dTT]) and distribution func- 
tions f±s{x,E) should be found from the kinetic equations, which follow from the 
Keldysh component of Eqs. ([l]) and for E >0 have a simple form d^f±s = within 
our approximations. These equations have no bound solutions: both distribution 
functions f±s{x,E) grow linearly with x far from the junction. Such a growth is 
limited in practice by inelastic collisions, which provide the spatial relaxation of 
f±s{x,E) to the equilibrium values at x ^ l± ^ ^o, where l± = \/Dz± are the 
inelastic scattering lengths and T± are the inelastic scattering times. To simplify 
the problem, instead of including complicated inelastic collision integrals, we add 
collision terms in the relaxation time approximation to the kinetic equations, 

lld^f+S = {f+S'feq)N{E), (13) 

fd^f^s=f~s/N{E), (14) 

where N{E) = Re{us) is the BCS DOS. 

We should supplement Eqs. ( I13b -( ll4b with proper boundary conditions on both 
left and right interfaces. On the tunnel barrier (x = 0) they follow from the Keldysh 
component of Eq. (O, 

dxf+s 1^=0 = (f+so -f+N): (15) 

1-0 = ^;^MiE) - ' ^^^^ 

where f±so{E) are the boundary values of f±s{x,E) at jc = 0. On the right trans- 
parent interface distribution functions become equilibrium functions of the right S 
reservoir, 

f+s\^^l^= feq, /-s|j.=^=0. (17) 

3 Single particle current 

The equation for the electric current follows from Eqs. O, ( 116b and for single 
particle current reads, 

I=^I^N{E){f.N-f-so)dE. (18) 
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Fig. 2 (Color online) IV characteristic of the NIS'S junction. = Ts = 03 Tc,W =0.01. Solid 
black line is calculated by use of Eq. i22\ . Other lines are calculated by use of Eq. <21t . L = 
50£,o: I- = 10^0 (red line); /_ = 20^o (green line); /_ = 50i^o (blue line); /_ = lOOi^o (dashed 
black line). 



To obtain /_5o we should solve the boundary problem (fT4]|.(fT6]l and (fTTl i. Doing 
this we get the following result, 

CC— 

^ a- + VmcoMp./^)) ' ^'^^ 
a-=2Wl-/^o=R{l-)/R, I5-=L/1-, (20) 

where R{1-) = l~/gN is the normal resistance of the superconducting lead per 
length From Eq. ( 1191 ) we see that the nonequilibrium correction to the current 
is of the order of a_ ^ W, which justifies our assumption about neglecting terms 
of the order of W in the kinetic equation. 

Substituting Eq. (fT9] l into Eq. dlSl l we finally obtain, 



/ = -L rNiE)f., c^iPjvm),E. (21) 

eRJA a_/v^V(£)+coth(/3_/v^VCE)) 

For L = or /_ = the boundary value f-so in Eq. ( flSl) is equal to zero and 
Eq. (1211) reduces to the well-known equation for single particle current between N 
and S reservoirs, connected through an insulating layer, 

/=— / N{E-eV)[nFiE-eV)-nF{E)]dE, (22) 

where nF{E) = [1 +exp(ii/r)]~^ is aFermi function. 

In the absence of inelastic relaxation in the S' lead (j3_ ^ 1) we can approxi- 
mate Eq. ^T\\ by 

l=^rf..^;^dE, (23) 
cRJa N{E) + aL 
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where = 2WL/^q = R{L) /R and R{L) = L/g^ is the resistance of the S' lead 
in the normal state. At zero temperature at eV ^ A the current given by Eq. ( 1231) 
can be calculated to first order in the small parameter A/eV I as follows, 

1 N^(E)dE A fV/^ x^dx V 
J = / i i = / ^ PS I / (24) 

eRjA N{E) + aL eRJi xV^^ + a^x^ ~ I) R,ot 

where = R + R{L) is the net normal resistance of the junction and the S' lead 
and lexc is an excess current given by the relation, 



eR 



ai 2af 

arctan 



1-«L (l-a2)'/' 




(25) 



Thus the IV characteristic Eq. ( 1231) exhibits a voltage-independent excess current 
at large voltage, which is the manifestation of the nonequilibrium in the S' lead. 
Here we should mention that the total excess current measured in the experiment is 
known to consist of the two contributions: the one coming from the single particle 
current at large voltage just calculated above, and the other coming from the two 
particle current (Andreev current). In this paper we do not calculate the latter 
contribution since it is of the order of W ^ a^. 

In Fig.|2]we plot the IV characteristic of the NIS'S junction for different values 
of I- parameter. I(V) given by Eqs. ( 121b . ( 1221) is an odd function of voltage and we 
plot it only for positive voltages. We fix 7/v = = 0.37^, where 7c is the critical 
temperature of the superconductor. For Aluminum 7]v = ~ 360 mK. We see that 
with the growth of the charge imbalance relaxation length Z_ the electric current 
decreases. When Z_ > L the length of the S' lead L plays the role of a characteristic 
relaxation length and the current is almost independent of Z_ . 



4 Cooling power 

The equation for the cooling power follows from Eqs. (H), ( fTSl) and reads 

P=-'V-^ ^f+N - f+so) dE, (26) 

where / is given by Eq. ( 1211) . To obtain f+so we should solve boundary problem 
(fT3]).(fT5]) and dTT] ). Doing this we get the following result. 



r ^ f+N a+ VNjE) + coth (/3+ y/NjE)) .27) 

a+^^V(£)+coth(j3+v^) ' 
a+=2Wl+/^o=R{l+)/Ry^W, ^+=L/l+, (28) 

where /?(/+) = l+/gN is the normal resistance of the superconducting lead per 
length /+. 

Substituting Eq. ( l27l) into Eq. ( 1261 ) we finally obtain, 

e^RjA ^ "^■'^^ ^ ^ a+V^VCE)+coth(/3+V^VCE)) 
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Fig. 3 (Color online) P(V) dependence of the NIS'S junction. TV = Ts = 0.3 Tc, W = 0.01. 
Solid black line is calculated by use of Eq. l |3Qt . Other lines are calculated by use of Eq. ( |29t . 
L = 50(^0 : /-=/+ = lOi^o (red line); /_=/+ = 20£,o (green line); L = /+ = 50i^o (blue line); 
L = /+ = 100^0 (dashed black line). 



For L = or Z-t = the boundary value /+50 in Eq. ( [26] l is equal to /^^^ and Eq 
reduces to the well-known equation for the heat current between N and S reser- 
voirs, connected through an insulating layer [6], 

1 /■+~ 

P=-^ N{E){E-eV)[nf{E-eV)-nF{E)]dE. (30) 

In the absence of inelastic relaxation in the S' lead (jS-t ^ 1) we can approxi- 
mate Eq. (|29ll by 

^--^^-^ .L ^^"-^ ' ^''^ l + a,N{Ey ^^^^ 

where the current / is given by Eq. ( |23l l. This equation corresponds to the case 
when the length of the S' lead L is smaller than the inelastic relaxation length and 
all relaxation occurs only in the S reservoir. 

In Figs. |3|4|5| we plot the P(V) dependence of the NIS'S junction for different 
values of l± parameters. P(V) given by Eqs. ( 1291 ), ( 130b is an even function of volt- 
age and we plot it only for positive voltages. From Fig. [3] it can be seen that with 
the growth of the relaxation lengths l± the cooling power decreases. Also the max- 
imal cooling power shifts to the region of smaller voltages than in the equilibrium 
case. We ascribe this suppression to the effect of backscattering on impurities and 
tunneling of nonequilibrium quasiparticles back to the normal metal reservoir. We 
can see that for large values of l± the cooling power is negative for all voltages. 

In Fig.|4]we plot P(V) for different ratio For a fixed Z+ length we vary 

the charge imbalance relaxation length /_. It can be seen that the cooling power 
increases with the growth of Z_ . This happens because of the decrease of the term 
IV in Eq. ( 1291) due to the suppression of the electric current (see Sec. [3]). 
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Fig. 4 (Color online) P(V) dependence of the NIS'S junction. Tn = Ts = 0.3 Tc, W = 0.01. Solid 
black line is calculated by use of Eq. i30\ . Other lines are calculated by use of Eq. <29t . L = 
50^0, /+ = W^o- 1- = lo|o (red line); L = 20£,o (green line); L = 50(^o (blue line); L = lOOi^o 
(dashed black line). 




Fig. 5 (Color online) P(V) dependence of the NIS'S junction. = Ts = 0.3 Tc, W = 0.001. 
Solid black line is calculated by use of Eq. l [30t . Other lines are calculated by use of Eq. l |29t . 
L = 50<^o: /-=/+ = lOi^o (red line); /_=/+ = 20^o (green line); /_=/+ = 50i^o (blue line); 
/_=/+ = 100(^0 (dashed black line). 



Finally we want to stress here the role of the transparency parameter W. In 
Fig. |5] we plot the P(V) dependence for a different value of tunneling parameter 
W = 10-3, which is one order of magnitude smaller than W in Figs. l3l4l Here we 
again see the decrease of the cooling power with the growth of l±, but the effect is 
smaller than in Fig. [3] This is obvious since the amplitudes of the nonequilibrium 
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Fig. 6 (Color online) Distribution functions f^s{E,x) (left plot) and Sf+s{E,x) = feq — f+s 
(right plot). Here Tn = Ts = 0.3 Tc, eV = 0.84, L = 50|o, l-=l+ = 10.^o, = 0.01. 

distribution functions f±sQ scale with the W parameter. For very strong tunnel 
barriers the nonequilibrium effect is therefore negligible. 

In order to estimate the characteristic parameters of the junctions for the trans- 
parency parameters W = 10^^ and W = 10^^ used in our numerical calculations, 
we will assume the junction area to be 200 x 200 nm and the thickness of the 
leads as well as the mean free path to be 50 nm. For Al leads, this results in the 
sheet resistance /?□ w 0.3 Q. and R{£,q) w 0.45 12 at ~ 300 nm. Then, accord- 
ing to Eq. (in, the tunneling probability and the junction resistance approach the 
values r w 2 X 10-^ R w 22.5 QforW = lO'^, and F w 2 x lO"'*, /? w 225 12 
for W = 10~^, respectively. 



5 Nonequilibrium quasiparticle distribution in S' 

Solving Eqs. (fT4]|.(fT6]l and (fTTl i we can obtain the function f-s{x,E) in the S' lead. 
It reads 



This equation describes the relaxation of the imbalance between the electrons and 
holes in the S' lead (branch mixing). From Eq. ( 1321) we see the exponential decay 
of the imbalance function /_ at the distance x ^ /_ y^N{E) from the tunnel barrier. 
AtE A the decay length diverges. 



Similarly from Eqs. (113b . (115b and ( 117b we obtain the function f+s{x,E) in the 
S' lead. 



f-NCC. 



sinh {p^-x/l_)/y^N{E) 



f-s 




(32) 



{f,g-f+N)a+^/N(E) 



sinh {l5+-x/l+)^/NiE) 



f+s — fe 



eq 




^. (33) 
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From this equation we see the exponential decay of the 5f+s function at the dis- 
tance X ^ Ij^I \JN{E) from the tunnel barrier, where 5f+s = feq — f+s- 

We note that the charge imbalance function /_ can be measured in the ex- 
periment [23]. In the adopted geometry. Fig. [H one can probe the /_s function 
at different locations in S' by small normal metal electrodes Nj, connected to S' 
through the insulating barrier at different points x,. Importantly, the area of these 
probing N'lS' junctions should be much smaller than the area of the initial NIS' 
barrier. Driving the current through the probing junction one can measure the cor- 
responding IV curve and calculate the function fs in the S' lead. 

In Fig. |6]we plot both f-s{E,x) and 5f+s{E,x) functions. We can see that 
nonequilibrium distributions are nonzero only at energies rather close to A . For 
large energies i? ^ Zi both functions are approaching zero. 

In this paper we used simple relaxation approximation approach. At zero tem- 
perature it gives qualitatively correct estimation of nonequilibrium properties of 
considered systems (see, for example, Ref. [24]). We note that this approximation 
leads to an incorrect expression of the charge imbalance relaxation length near the 
critical temperature Ii25i] . At intermediate temperatures Q <T <Tc this approxi- 
mation gives a result correct with an accuracy of a numerical factor. Since in this 
paper we consider the temperatures = Tg equal to 0.37c which is much smaller 
than Tc our description is qualitatively correct. However, we should stress that 
for higher temperatures one should add inelastic collision integrals in the kinetic 
equations [HI, 111 HI- 

6 Conclusion 

We investigated the electric and heat currents in an NIS'S junction in the diffusive 
limit. We have developed a model which describes the nonequilibrium quasiparti- 
cle injection and relaxation in the superconducting lead. This model will be used 
as a tool to fit experimental data in various types of NIS tunnel junctions in non- 
reservoir geometries. These fits will be presented elsewhere. 

We showed that in the case when relaxation lengths in the superconductor are 
rather long compared to the coherence length, the electric current and the cooling 
power for eV below A are suppressed. We ascribe this suppression to the backtun- 
neling of nonequilibrium quasiparticles into the normal metal. The value of this 
suppression scales with the diffusive transparency parameter W. 

Finally, we calculated the nonequilibrium distribution functions in the super- 
conducting lead. 
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